Permeabilities of metamaterials 
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Abstract 



Scattering of electromagnetic (EM) waves by many smaU particles, 
Q^. embedded in a given medium, is studied. Physical properties of the 

Q I particles are described by their boundary impedances. The limiting 

. equation is obtained for the effective EM field in the limiting medium, 

O ' in the limit a — >■ 0, where a is the characteristic size of a particle and 

c/3 i the number AI{a) of the particles tends to infinity at a suitable rate. 

An analytical formula for the permeability ^(a;) of the limiting medium 
is given. Analysis of this formula allows one to find out the range of 
the values of the permeability in the material, obtained by embedding 
many small particles. 
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. 1 Introduction 

In this paper we outline a theory of monochromatic electromagnetic (EM) 
^ , wave scattering by many small particles (bodies) embedded in a homoge- 

^ I neous medium which is described by the permittivity e' = eo + i^' > 0, 

and constant permeability hq. The wave number k = a;(e'/Lio)^^^- Related 
papers are [l]-[8], and the results we use without proof are taken mostly 
from 0]. The small particles are embedded in a finite domain fi. Smallness 
of the particles means that \k\a <ti 1, where a is the characteristic size of 
the particles, and k is the wave number. The medium, created by the em- 
bedding of the small particles, has new physical properties. In particular, 
it has a spatially inhomogeneous magnetic permeability n{x), which can be 
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controlled by the choice of the boundary impedances of the embedded small 
particles and their distribution density. This is a new physical effect. An 
analytic formula for the permeability of the new medium is derived in |4j: 

/^(-) = (1) 

where 

^^x) = l + l^h{x)N{x). (2) 

Here u is the frequency of the EM field, eo is the constant dielectric pa- 
rameter of the original medium, h(x) is a function describing boundary 
impedances of the small embedded particles, and N{x) > is a function 
describing the distribution of these particles. We assume that in any sub- 
domain A, the number Af{A) of the embedded particles Dm is given by the 
formula: 

A/'(A) = / N{x)dx[l + o{l)], a^O, (3) 

where N{x) > is a continuous function, vanishing outside of the finite 
domain 0, in which small particles (bodies) are distributed, 1 < m < M, 
M = M{a), K € (0, 1) is a parameter, and the boundary impedances of the 
small particles are defined by the formula 

_ h{Xm) p r, (A) 

where point inside m— th particle Dm, Re h{x) > 0, and h{x) is a 

continuous function vanishing outside fi. The impedance boundary condition 
on the surface Sm of the m— th particle Dm is = C,m[H^-,N]-, where £"* 
{H^) is the tangential component of E (H) on Sm, and N is the unit normal 
to Sm, pointing out of Dm- Physical properties of the impedance C are 
discussed in [2]. In particular, Re^ > 0. 

Since one can choose the functions N{x) and h{x), one can create a 
desired magnetic permeability in 0. This is a novel idea, to the author's 
knowledge. 

We have also derived in [3] an analytic formula for the refraction coeffi- 
cient of the medium in created by the embedding of many small particles. 

KHx)= ,,,, , , k^ = u'e'fio, (5) 

where the coefficient appears if Dm are balls of radius a centered at the 
points Xm- For the small bodies Dm of arbitrary shape the coefficient % 
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should be replaced by a tensorial coefficient Cm = 0{1), depending of the 
shape of Dm- 

An equation for the EM field in the limiting medium is derived in [3]: 

E{x) = Eo{x) - -^V X / g{x,y)h{y)N{y)V x E{y)dy. (6) 

Jn 

The limiting medium is created when the size a of small particles tends to 
zero while the total number M = M{a) of the particles tends to infinity at 
the rate determined by ([3]). 

The refraction coefficient in the limiting medium is spatially inhomoge- 
neous. 

Our theory may be viewed as a " homogenization theory", but it dif- 
fers from the usual homogenization theory (see, e.g., [l], [3], and references 
therein) in several respects: we do not assume any periodic structure in the 
distribution of small bodies, our operators are non-selfadjoint, the spectrum 
of these operators is not discrete, etc. Our ideas, methods, and techiques 
are quite different from the usual methods. 

These ideas are similar to the ideas developed in papers [Sj [6] , where 
scalar wave scattering by small bodies was studied, and in the papers [3] 
and [8], where EM wave scattering by many small particles was stidied. 
Scattering of EM waves brought new technical difficulties which were re- 
solved in [4] and [8]. The difficulties come from the vectorial nature of the 
boundary conditions. Results from [7] were used for a justification of the 
passages to the limit a ^ 0. 

In [1] a new numerical method for solving many-body wave-scattering 
problems for small scatterers is given. 

We predict theoretically the new physical phenomenon: by embedding 
many small particles with suitable boundary impedances into a given homo- 
geneous medium, one can create a medium with a desired spatially inhomo- 
geneous permeability ([1]). 

One can create material with a desired permeability ^i{x) by embedding 
small particles with suitably chosen boundary impedances. Indeed, by for- 
mula ([T]) one can choose a complex- valued, in general, function /i(x), and a 
non-negative function N{x) > 0, describing the density distribution of the 
small particles, so that the right-hand side of formula ([T]) will yield a desired 
function fi{x). 

The goal of this paper is to discuss the possible values of magnetic per- 
meability one can get in the metamaterial, obtained by embedding many 
small impedance particles in a given material. 

This discussion is given in the next Section. 
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2 Possible values of permeability 



Let us denote q = q{x) := and rewrite formula ([T|) as 



/io(l - qh2 - iqhi) 
(1 - Q'n2)'' + q^hl 



Quantity q{x) > 0. Thus, the values of n{x) can be controlled by choosing 
the function h in equation 
Conclusion: 

Assume that the impedances of the small particles are defined by formula 
dll, where h[x) = hi{x) + ih2{x) is a function with Reh[x) = hi[x) > 0. 
Let fi{x) = fJ,i{x) + ifi2{x), where ^i(x) is the real part of fi{x). Then it 
is not possible to create metamaterial with f^2ix) > by embedding many 
small impedance particles into a given material. Indeed, 

(1 - qh2{x)Y + q^hf{x) 

provided that /xq, q and hi are all positive, which is the case. The assumption 
hi{x) > is a necessary assumption for boundary impedances, see f^, p. 301. 

For high frequencies u the quantity q{x) = 0{^) is small, so if hi{x) = 
0(1), then fj.2{x) = O(^) is also small. 
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